Calibration of frequency monitors having dual etalon signals in quadrature

ABSTRACT

Improved calibration of a dual-etalon frequency monitor having x-y outputs is provided. An ellipse is fit to the (x,y) points from a set of calibration data. For each (x,y) point, an angle θ is determined. A linear fit of frequency to θ is provided. Differences between this linear fit and the determined values of θ are accounted for by including a spline fit to this difference in the calibration.

CROSS REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. provisional patent application No. 60/873,652, filed on Dec. 7, 2006, entitled “Laser wavelength targeting and auto-recalibration of a dual-etalon wavelength monitor”, and hereby incorporated by reference in its entirety.

FIELD OF THE INVENTION

This invention relates to calibration of optical frequency or wavelength monitors.

BACKGROUND

Measurement of frequency (or wavelength) of an optical source is generally useful in a wide variety of situations, and is particularly important in applications where very precise control of laser output frequency is required, such as high precision spectroscopy. Accordingly, various methods for optical frequency monitoring have been developed to date. One such method is the use of one or more etalons to provide suitable signals for frequency monitoring. Etalons are attractive for frequency monitoring because they can provide high accuracy and precision without being unduly large or expensive.

One type of etalon frequency monitor relies on the use of two etalon signals in quadrature, because the use of two quadrature signals helps avoid reduced measurement sensitivity at local maxima or minima of the etalon signals. One example of this kind of approach is described in U.S. Pat. No. 6,859,284 and in U.S. Pat. No. 6,331,892, where detectors are disposed in an etalon interference pattern at a quadrature separation from each other. In U.S. Pat. No. 6,178,002, two beams take slightly different paths through an etalon, where the difference in path length is selected such that signals corresponding to the two beams are in quadrature.

The specific issue of how to derive a frequency reading from measured quadrature etalon signals has also been investigated. For example, in US 2007/0195328, “dead zones” of reduced sensitivity for each etalon signal are identified. If an etalon signal is in its dead zone, it is not included in the frequency computation. In this manner, the frequency computation is based on one or both of the etalon signals, and relatively useless etalon signals are systematically excluded from the computations, thereby improving performance.

However, the approach of defining dead zones in this manner can encounter significant difficulties in practice. In particular, switching from the use of one etalon signal to the use of two etalon signals at a dead zone boundary can lead to an undesirable situation where the frequency monitor error (i.e., the difference between monitor readout and actual frequency) is discontinuous. Although such discontinuities can be removed in calibration, they can reappear as a result of calibration drift. A discontinuous monitor error is particularly undesirable in situations where the monitor is part of a frequency control loop, since the loop may become unstable, or it may be impossible to access the target frequency at all.

Accordingly, it would be an advance in the art to derive frequency readings from a frequency monitor having two etalon signals in quadrature in a manner that provides high precision and accuracy and also inherently provides a continuous monitor error as a function of frequency.

SUMMARY

Improved calibration of a dual-etalon frequency monitor having x-y outputs is provided. An ellipse is fit to the (x,y) points from a set of calibration data. For each (x,y) point, an angle θ is determined. A linear fit of frequency to θ is provided. Differences between the calibration data and the linear fit are accommodated by including a spline fit to this difference in the calibration.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 shows an example of two measured etalon signals in quadrature from a frequency monitor.

FIG. 2 shows an example of an elliptical fit to etalon signals as employed in an embodiment of the invention.

FIG. 3 shows an example of the relationship between frequency and angle variable and its proximity to a linear function as employed in an embodiment of the invention.

FIG. 4 shows the difference between the calibration data and the linear fit of the example of FIGS. 1-3.

FIG. 5 shows a spline fit to the difference shown in FIG. 4 as employed in an embodiment of the invention, over a reduced angular range.

DETAILED DESCRIPTION

In order to better appreciate the present invention, it is helpful to consider how quadrature etalon signals have previously been employed to compute wavelength. In U.S. Pat. No. 6,331,892, it is assumed that the two etalon signals S1 and S2 are normalized to have equal amplitude and are proportional to sin(2πλ/Λ) and cos(2πλ/Λ), where λ is wavelength and Λ is the period of the etalon signals. Under these conditions, tan⁻¹(S1/S2) and cot⁻¹(S2/S1) are both proportional to wavelength, so with an appropriate linear scaling these quantities can provide the desired wavelength readout. A similar approach is considered in U.S. Pat. No. 6,178,002. As is well known in the art, etalon signals tend to be approximately periodic in both frequency and wavelength, so this kind of approach would also be applicable to measurement of frequency.

However, this relatively simple computational approach can lead to a significant loss of accuracy in the computed wavelength (or frequency) readout in practice. According to embodiments of the present invention, improved accuracy is provided by a calibration method including the following steps: a) providing calibration data; b) fitting an ellipse to the calibration data points; c) determining an angle parameter of each calibration data point from the fitted ellipse; d) determining a linear fit of frequency (or wavelength) to the angle parameter; e) determining the difference between the calibration data and the linear fit; f) fitting a spline to the difference; and g) providing the ellipse fit, the line fit and the spline fit as calibration outputs. These steps are explained in greater detail by way of the following example.

FIG. 1 shows an example of measured frequency monitor data. Here laser frequency is measured by a wave meter, and the WLM ratio signals are outputs from a typical dual-etalon frequency monitor. It is convenient to regard these calibration data as being in the form of triples (x_(i), y_(i), ν_(i)) indexed by i. For each triple, ν_(i) is a known frequency, x_(i) is a first etalon signal at ν_(i), and y_(i) is a second etalon signal at ν_(i). The first and second etalon signals are substantially in quadrature (i.e., have a phase difference of about π/2, as shown in FIG. 1). Preferably, the etalon signals x and y are ratio signals that are substantially independent of the optical power received by the frequency monitor. Suitable arrangements of optical components for providing such etalon signals are well known in the art. For brevity it is convenient to regard any frequency monitor providing two etalon signals in quadrature as being a “dual-etalon frequency monitor”, even if the monitor includes one etalon as opposed to two distinct physical etalons.

FIG. 2 shows a plot of the (x_(i), y_(i)) points in the x-y plane. In this example, the points tend to lie along an ellipse. For ideal etalon signals, the (x_(i), y_(i)) points should fall exactly on an ellipse. It is convenient to describe this ideal ellipse with the following equations:

x=x ₀ +A cos θ, and   (1)

y=y ₀ +B sin(θ+ε).   (2)

Here parameters x₀ and y₀ are the center of the ellipse, and parameters A and B account for the possibility of different amplitudes of x and y. The parameter ε accounts for the possibility of imperfect quadrature (i.e., x and y having a phase difference not exactly equal to an odd multiple of π/2). Ideally, the angle parameter θ is exactly proportional to frequency. Thus calibrating a dual-etalon frequency monitor can be regarded as determining the best linear relation between θ and ν, as well as providing corrections to this linear relation. Note that θ is not the geometrical ellipse angle unless ε=0. Allowing for unequal x and y amplitude and imperfect quadrature as in Eqs. 1 and 2 tends to provide substantially improved accuracy compared to the simpler approach of U.S. Pat. No. 6,331,892, which corresponds to assuming A=B and ε=0 in Eqs. 1 and 2.

Although any method of fitting an ellipse to the (x_(i), y_(i)) points can be employed in practicing embodiments of the invention, it is preferred to begin by fitting the (x_(i), y_(i)) points to a quadratic form

a ₀ x ² +a ₁ xy+a ₂ y ² +a ₃ x+a ₄ y+a ₅=0,   (3)

to determine fitting parameters a₀, a₁, a₂, a₃, a₄, and a₅. Methods for performing such fitting are known, e.g., as described by Fitzgibbon et al. in “Direct Least Squares Fitting of Ellipses”, IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 21, no. 5, pp. 476-480, May 1999. Once the parameters a₀, a₁, a₂, a₃, a₄, and a₅ are known, the parameters of interest (i.e., x₀, y₀, A, B, and ε) can be determined as follows.

First, x₀ and y₀ are given by

$\begin{matrix} {x_{0} = {{\frac{{a_{1}a_{4}} - {2a_{2}a_{3}}}{{4a_{0}a_{2}} - a_{1}^{2}}\mspace{14mu} {and}\mspace{14mu} y_{0}} = {\frac{{a_{1}a_{3}} - {2a_{0}a_{4}}}{{4a_{0}a_{2}} - a_{1}^{2}}.}}} & (4) \end{matrix}$

Next, it is convenient to define r₁ and r₂ to be the lengths of the principal semi-axes, and φ as the angle which the first principal axis makes with respect to the x axis. These quantities are given by

$\begin{matrix} {{{\tan \; 2\varphi} = \frac{a_{1}}{a_{0} - a_{2}}}{r_{1} = \sqrt{\frac{{a_{0}x_{0}^{2}} + {a_{1}x_{0}y_{0}} + {a_{2}y_{0}^{2}} - a_{5}}{{a_{0}\cos^{2}\varphi} + {a_{1}\cos \; {\varphi sin\varphi}} + {a_{2}\sin^{2}\varphi}}}}r_{2} = \sqrt{\frac{{a_{0}x_{0}^{2}} + {a_{1}x_{0}y_{0}} + {a_{2}y_{0}^{2}} - a_{5}}{{a_{0}\sin^{2}\varphi} - {a_{1}\cos \; \varphi \; \sin \; \varphi} + {a_{2}\cos^{2}\varphi}}}} & (5) \\ {{Finally},{{we}\mspace{14mu} {have}}} & \; \\ {{A = \sqrt{{r_{1}^{2}\cos^{2}\varphi} + {r_{2}^{2}\sin^{2}\varphi}}}{B = \sqrt{{r_{1}^{2}\sin^{2}\varphi} + {r_{2}^{2}\cos^{2}\varphi}}}{{\cos \; ɛ} = \frac{r_{1}r_{2}}{AB}}} & (6) \end{matrix}$

where the sign of E is negative if (r₁−r₂)sin φ cos φ<0. Following this procedure for the example of FIG. 2 gives x₀=0.8316, y₀=0.8246, A=0.8211, B=0.8001, and ε=−0.0513. The resulting ellipse is shown with a solid line in FIG. 2.

Once the best fit ellipse is known, an angle parameter θ_(i) can be associated with each (x_(i), y_(i)) point. Since the (x_(i), y_(i)) points do not necessarily lie on the best fit ellipse, it is usually not possible to find a θ_(i) that satisfies Eqs. 1 and 2 simultaneously. There are various reasonable ways to calculate θ_(i) from (x_(i), y_(i)), any of which can be employed in practicing the invention. Preferably θ_(i) is determined from

θ_(i) =a tan 2(A(y _(i) −y ₀)−B(x_(i) −x ₀)sin ε, B(x _(i) −x ₀)cos ε),   (7)

where a tan 2(y,x) is the arctangent of y/x with a range of (−π,π] (i.e., a four quadrant angle). This procedure is equivalent to drawing a line from the center of the best fit ellipse to (x_(i), y_(i)) and computing θ_(i) using Eqs. 1 and 2 as applied to the intersection point of this radial line and the ellipse.

The result of Eq. 7 is an angle parameter in the range (−π,π], but the frequency range of interest may cover a larger range of angles. In other words, there may be multiple revolutions along the ellipse, e.g., as shown in FIGS. 1 and 2. In such cases, it is convenient to add multiples of 2π as necessary to θ_(i) as determined by Eq. 7 in order to establish a continuous relation between frequency and angle, e.g., as shown in FIG. 3. Since the absolute frequencies ν_(i) of the calibration data are known, this is a straightforward process.

As seen in FIG. 3, the relation between frequency and angle is typically highly linear. Accordingly, fitting a line ν=sθ+ν₀ to the set of points (θ_(i), ν_(i)) provides a linear angle-frequency fit, so the parameters s and ν₀ can be regarded as first order calibration parameters. Any method of linear fitting (e.g., linear least squares) can be employed to determine this linear fit.

However, for high precision applications, it can be important to account for differences between this linear fit and the calibration data. FIG. 4 shows a plot of this difference for the example of FIGS. 1-3. Here, the difference e_(i) of the linear angle-frequency fit is defined by e_(i)=ν_(i)−(sθ_(i)+ν₀). The relative scale of the vertical axes of FIGS. 3 and 4 is noteworthy. The maximum difference magnitude in FIG. 4 is about 0.015 cm⁻¹, while the frequency in FIG. 3 is on the order of 6550 cm⁻¹. Thus the differences shown in FIG. 4 are on the order of 1 part in 400,000 or less. High precision calibration to further reduce this already low difference can be crucial for high precision spectroscopy applications (e.g., cavity ring-down spectroscopy for trace gas detection), but is not required for many other applications (e.g., wavelength division multiplexing for telecommunication).

For calibration, it is helpful to use an appropriately smoothed fit to the differences e_(i) as the correction to the first order linear model. Smoothing helps to reduce the effect of random measurement noise on the calibration. Preferably, this smoothed representation is obtained by fitting a spline to the set of points (θ_(i), e_(i)) to provide a spline difference fit e(θ). More preferably, cubic B-splines having equally spaced knots are employed. For example, FIG. 5 shows the difference (circles) and corresponding B-spline fit (solid line) for part of the angular range of FIG. 4. In this example, the knots of the B-spline are spaced by π/8 radians. By fixing the knot spacing, the spline fit inherently performs a desirable smoothing function as compared to a spline constrained to pass through every point e_(i). The amount of smoothing can be adjusted by altering the predetermined knot spacing, since increasing the knot spacing increases the smoothing effect, and decreasing the knot spacing decreases the smoothing effect. Methods for performing such spline fits are well known in the art.

At the end of this process, a set of v(θ) calibration outputs including the ellipse parameters A, B, x₀, y₀, and ε, the line parameters s and ν₀, and the spline fit e(θ) is available for further use. The calibration curve is given by ν=sθ+v₀+e(θ).

The preceding description relates to a calibration method. Further aspects of embodiments of the invention relate to using the ν(θ) calibration outputs to measure the emission frequency of an optical source, to using the ν(θ) calibration outputs to set the emission frequency of an optical source, and to updating the ν(θ) calibration.

Once a v(θ) calibration is available, it can be employed to measure operating frequency of an optical source in a system including the calibrated dual-etalon frequency monitor. The measurement begins by measuring source operating data (x_(op), y_(op)), where x_(op) and y_(op) are the etalon signals provided by the calibrated frequency monitor responsive to radiation from the optical source. An operating angle parameter θ_(op) can be calculated according to

θ_(op) =a tan 2(A(y _(op) −y ₀)−B(x _(op) −x ₀)sin ε, B(x _(op) −x ₀)cos ε).   (8)

Note that the ellipse parameters of the calibration are required in order to determine θ_(op) from x_(op) and y_(op). From θ_(op), the operating frequency ν_(op) can be determined according to ν_(op)=sθ_(op)+ν₀+e(θ_(op)), and provided as an output.

In cases where the frequency range of interest corresponds to an angular range greater than 2π, it is necessary to remove the 2πn ambiguity. A relatively low precision frequency measurement will suffice to remove this ambiguity. More specifically, a set of angles θ_(n)=θ_(op)+2 πn can be defined, where n is any integer that results in θ_(n) being within the θ calibration range. A corresponding set of frequencies ν_(n)=ν(θ_(n)) can be calculated as described above. The frequency ν_(n) that is closest in value to the result of a secondary frequency measurement can be provided as the frequency output. Since the required precision for the secondary measurement is relatively low (i.e., it only needs to be able to distinguish ν_(n) from ν_(n±1)), it is not particularly difficult to provide such a secondary measurement, and any such method can be employed in practicing the invention. Note that e(θ) has a domain which includes the entire calibration range of θ, which can include many multiples of 2π.

Frequently, an indirect secondary measurement can be sufficiently accurate. In one system design, the etalon free spectral range is 1.6 cm⁻¹ (i.e., about 48 GHz), and the optical source is a distributed feedback (DFB) laser having a temperature tuning coefficient of 14 GHz/C. Successive revolutions of the ellipse are separated by about 3.4° C., so knowledge of the DFB temperature with about 1° C. precision or better can provide a sufficiently accurate secondary frequency measurement.

The v(θ) calibration as described above can also be employed for setting the operating frequency of an optical source to a desired value ν_(in) in a system having a calibrated dual-etalon frequency monitor. The first step here is to determine an angle θ_(in) corresponding to ν_(in) according to ν_(in)=sθ_(in)+ν₀+e (θ_(in)). Various methods can be employed for this calculation. One approach is inverse cubic interpolation to solve for θ_(in), which is appropriate because e(θ_(in)) is a third-degree function of θ_(in). Another approach is iterative numerical solution of ν_(in)=sθ_(in)+ν₀+e (θ_(in)), using θ*=(ν_(in)−v₀)/s as a starting point. In cases where the frequency range of interest corresponds to an angular range greater than 2π, coarse frequency control of the optical source sufficient to distinguish ν_(in) from any frequency ν_(n)=ν(θ_(in)+2πn), where n is any non-zero integer, can be provided. Such frequency control is typically not difficult, and any such method can be employed in practicing the invention.

For fine frequency control, it is necessary to control the optical source such that the etalon signals x and y correspond to the angle θ_(in). It is important to note that such control should not attempt to enforce Eqs. 1 and 2 simultaneously, since it is ordinarily not possible to satisfy both these equations simultaneously. Recall in FIG. 2 that most of the data points are not on the calibration ellipse defined by Eqs. 1 and 2. Instead a control condition of the form

w ₁(x−x ₀)+w ₂(y−y ₀)=0   (9)

can be enforced, where the weights w₁ and w₂ are preferably each in the range (−1, 1]. A preferred choice for these weights is given by

$\begin{matrix} {{w_{1} = {{- K}\; \frac{\sin \left( {\theta_{in} + ɛ} \right)}{A\; \cos \; ɛ}}}{w_{2} = {K\; \frac{\cos \; \theta_{in}}{B\; \cos \; ɛ}}}} & (10) \end{matrix}$

where K is an overall scale factor that can be employed to appropriately scale w₁ and w₂. This choice is preferred because it causes the derivative of the left hand side of Eq. 9 with respect to θ_(in) at the desired set point to be equal to −K, which is independent of θ_(in) and all the parameters of the wavelength monitor, which is highly desirable for the control loop. Since the ratio of w₁ to w₂ is the important quantity for setting θ, the scaling can be selected for convenience of control loop implementation. For a reasonable amount of calibration drift (i.e., a drift which is not so large as to cause the nominal ellipse center (x₀, y₀) to lie outside the actual path defined by the (x_(i), y_(i)) points), the locking condition of Eq. 9 can be satisfied.

In general, any set of known frequencies having corresponding measured etalon signals can be employed to generate an initial calibration and/or to update an existing calibration. Calibration update can be full or partial. A full calibration update would include re-calculation of the ellipse parameters based on the new calibration data, while a partial calibration update would make use of the existing ellipse parameters and update the linear fit and/or the spline fit.

In many cases, it is preferred to perform partial calibration updates incrementally, e.g., by making incremental changes to the coefficients of the spline fit e(θ). Since calibration updates typically result in small changes to the spline fit coefficients, computation time can be significantly reduced by using the existing spline coefficients as a starting point in an iterative update calculation. Methods for such iterative updating of a spline fit are well known in the art, and any such method can be employed in practicing the invention.

An output frequency error Δν_(j)=ν_(j,out)−ν_(j) can be measured at several points, where ν_(j) is the frequency according to the calibrated system (i.e., related to θ according to the ν(θ) calibration), and ν_(j,out) is the corresponding measurement of the source frequency (e.g., as provided by an independent instrument, such as a wave meter). The spline difference fit e(θ) can be revised to decrease the output frequency error (e.g., by reducing the squared error). Preferably, such revision is performed according to an iterative procedure to minimize the output frequency error.

In some applications of high precision frequency monitoring, there are special features of the overall system that can be exploited to expedite calibration. For example, in a cavity ring-down spectroscopy (CRDS) instrument, the ring-down cavity can serve as a calibration frequency measurement standard because the resonant frequencies of the cavity are equally spaced by the free spectral range (FSR) of the cavity. There are several ways this feature can be exploited for calibration.

For example, the frequencies for initial calibration or calibration update of the frequency monitor can be selected to be equally spaced by the ring-down cavity FSR. This possibility is particularly convenient for calibration update, because it is possible to perform ring-down measurements in such a way that helpful data for calibration update are automatically generated. In a typical CRDS instrument, the length of the ring-down cavity can be varied over a short distance (i.e., on the order of an optical wavelength) by a suitable mechanical actuator, such as a piezoelectric transducer (PZT). The PZT is typically scanned repetitively through a fixed range of motion in order to generate ring-down events when a mode of the ring-down cavity coincides in frequency with the source frequency.

It is possible to adjust the nominal frequencies of a source operating frequency scan such that the PZT voltages are equal for all ring-down events. For example, the nominal separation between frequencies of the scan can be adjusted to remove any systematic dependence of PZT voltage on frequency, and then the scan frequencies can be individually adjusted as needed to equalize the PZT voltages as much as possible.

When this condition of equal PZT voltages holds, the frequency spacing for the source tuning is equal to the ring-down cavity FSR, because the cavity length is the same at all ring-down events. In this situation, the calibration can be updated because the actual frequencies of the ring-down events are known to be exactly equally separated in frequency. A further advantage of this arrangement is that the resulting calibration update does not depend at all on the PZT sensitivity (i.e., the μm/V provided by the PZT), since all data are obtained at the same PZT voltage.

Such information can be employed to perform a relative calibration update or an absolute calibration update. An absolute calibration update can be performed if the cavity FSR is known and if the frequency of one of the calibration update points is independently known (e.g., by reference to a spectral line). Such absolute calibration updating can proceed as described above.

For a relative calibration update, it is preferred to impose a constraint of no average frequency shift in the update. More specifically, suppose an equal PZT voltage FSR measurement scheme as described above provides a set of (x_(i), y_(i), θ_(i), ν_(i)) points, where x_(i) and y_(i) are the etalon signals, θ_(i) is the angle corresponding to x_(i) and y_(i), and ν_(i) is the nominal frequency calculated using θ_(i) according to the existing v(θ) calibration data. Calibration drift can lead to a situation in which the ν_(i) are not equally spaced. The v(θ) calibration data can be updated such that a set of nominal frequencies ν_(i)′ corresponding to θ_(i) as calculated from the updated calibration data are equally spaced. As indicated above, it is preferred to constrain the updating method such that the average of the ν_(i)′ is equal to the average of the ν_(i). Methods for updating a spline fit subject to such a constraint are known in the art. The cavity FSR can be provided as an input to the updating method, but this is not essential, because the frequency spacing can also be regarded as a parameter of the updating fit.

The preceding description of relative calibration updates is based on setting the frequency spacing equal to the cavity FSR. In some cases where data points having different PZT voltages are employed, sufficient accuracy for calibration updating can be provided by including an estimate of the PZT sensitivity in the calculations.

The preceding description has been by way of example as opposed to limitation. In particular, the same ideas are applicable to calibration of a wavelength monitor, setting the operating wavelength of an optical source, updating wavelength monitor calibration, and measuring wavelength with a calibrated dual-etalon wavelength monitor. The basis for this equivalence is the approximate periodicity of etalon signals with respect to both frequency and wavelength. The signal from a non-dispersive etalon is periodic in frequency, so it can only be approximately periodic in wavelength because wavelength is inversely proportional to frequency. However, in many wavelength monitor applications, this lack of exact periodicity is not troublesome. 

1. A method of calibrating a frequency monitor having dual etalon signals in quadrature, the method comprising: a) providing calibration data in the form of triples (x_(i), y_(i), ν_(i)) indexed by i, wherein for each triple, x_(i) is a first etalon signal at a known frequency ν_(i) and y_(i) is a second etalon signal at said known frequency ν_(i), and wherein said first and second etalon signals are substantially in quadrature; b) fitting an ellipse to the set of points (x_(i), y_(i)), wherein said ellipse is parametrically defined according to x=x ₀ +A cos θ y=y ₀ +B sin(θ+ε); c) calculating an angle parameter θ_(i) of said ellipse for each of said points (x_(i), y_(i)); d) fitting a line to the set of points (θ_(i), ν_(i)) to provide a linear angle-frequency fit, wherein said line is parametrically defined according to ν=sθ+ν₀; e) providing said parameters s and ν₀ and said parameters A, B, x₀, y₀, and ε as ν(θ) calibration outputs.
 2. The method of claim 1, further comprising: f) determining a difference e_(i) of said linear angle-frequency fit at each of said triples according to e _(i)=ν₁−(sθ _(i)+ν₀); g) fitting a spline to the set of points (θ_(i), e_(i)) to provide a spline difference fit e(θ); h) providing said spline difference fit e(θ) as part of said v(θ) calibration outputs.
 3. The method of claim 1, wherein said calculating an angle parameter θ_(i) of said ellipse for each of said points (x_(i), y_(i)) comprises setting θ_(i)=a tan 2(A(y_(i)−y₀)−B(x_(i)−x₀)sin ε, B(x _(i)−x₀)cos ε).
 4. The method of claim 1, wherein said frequencies ν_(i) of said calibration data are equally spaced.
 5. The method of claim 4, wherein said frequencies ν_(i) of said calibration data are determined by a reference optical cavity.
 6. The method of claim 1, wherein said fitting an ellipse comprises: fitting a quadratic form a₀x²+a₁xy+a₂y²+a₃x+a₄y+a₅=0 to said set of points (x_(i), y_(i)) to provide values for fitting parameters a₀, a₁, a₂, a₃, a₄, and a₅; calculating values for x₀, y₀, A, B, and ε, from said fitting parameters a₀, a₁, a₂, a₃, a₄, and a₅.
 7. The method of claim 2, wherein said fitting a spline comprises fitting cubic B-splines defined on a uniformly spaced set of knots.
 8. A method of measuring an operating frequency of an optical source in a system including a frequency monitor having dual etalon signals in quadrature, the method comprising: providing v(θ) calibration outputs according to the method of claim 1; measuring source operating data (x_(op), y_(op)), wherein x_(op) and y_(op) are respectively said first and second etalon signals provided by said frequency monitor responsive to radiation from said optical source having said operating frequency; calculating an angle parameter θ_(op) of said optical source according to θ_(op)=a tan 2(A(y_(op)−y₀)−B(x_(op)−x₀)sin ε, B(x_(op)−x₀)cos ε); determining said operating frequency from said ν(θ) calibration outputs and from said angle parameter θ_(op); providing said determined operating frequency as an output.
 9. The method of claim 8, wherein said determining said operating frequency from said ν(θ) calibration outputs and from said angle parameter θ_(op) comprises: providing a secondary measurement of said operating frequency; defining a set of angles θ_(n)=θ_(op)+2 πn, where n is any integer that results in θ_(n) being within a predetermined calibration range; calculating a set of frequencies ν_(n)=ν(θ_(n)) from said ν(θ) calibration outputs and from said set of angles θ_(n); selecting one of said set of frequencies ν_(n) as said determined operating frequency according to said secondary measurement.
 10. A method of setting an optical source to operate at or near a predetermined frequency in a system including a frequency monitor having dual etalon signals in quadrature, the method comprising: providing ν(θ) calibration outputs according to the method of claim 2; determining an angle θ_(in) corresponding to said predetermined frequency ν_(in) from said ν(θ) calibration outputs; controlling the output frequency of said optical source such that w₁(x−x₀)+w₂(y−y₀)=0, wherein x and y are said first and second etalon signals respectively, and wherein $w_{1} = {{{- \frac{\sin \; \left( {\theta_{in} + ɛ} \right)}{{Acos}\; ɛ}}\mspace{14mu} {and}\mspace{14mu} w_{2}} = {\frac{\cos \; \theta_{in}}{B\; \cos \; ɛ}.}}$
 11. The method of claim 10, further comprising: providing coarse frequency control of said optical source sufficient to distinguish said predetermined frequency ν_(in) from any frequency ν_(n)=ν(θ_(in)+2 πn), where n can be any non-zero integer; wherein said frequencies ν_(n) are determined from said corresponding angles θ_(in)+2 πn according to said ν(θ) calibration outputs.
 12. A method of calibration update in a system having a frequency monitor having dual etalon signals in quadrature, the method comprising: setting an optical source to operate at or near each of a set of predetermined frequencies ν_(j) according to the method of claim 10; measuring an output frequency error Δν_(j)=ν_(j,out)−ν_(j) at each of said set of predetermined frequencies, where ν_(j,out) is optical source frequency when said source is set to operate at ν_(j), wherein ν_(j,out) is measured independently of said ν(θ) calibration outputs; revising said spline difference fit e(θ) to decrease said output frequency error, by making incremental changes to spline fit coefficients.
 13. The method of claim 12, wherein said revising said spline difference fit is repeated one or more times.
 14. A method of calibrating a wavelength monitor having dual etalon signals in quadrature, the method comprising: a) providing calibration data in the form of triples (x_(i), y_(i), ν_(i)) indexed by i, wherein for each triple, x_(i) is a first etalon signal at a known wavelength λ_(i) and y_(i) is a second etalon signal at said known wavelength λ_(i), and wherein said first and second etalon signals are substantially in quadrature; b) fitting an ellipse to the set of points (x_(i), y_(i)), wherein said ellipse is parametrically defined according to x=x ₀ +A cos θ y=y ₀ +B sin(θ+ε); c) calculating an angle parameter θ_(i) of said ellipse for each of said points (x_(i), y_(i)); d) fitting a line to the set of points (θ_(i), λ_(i)) to provide a linear angle-wavelength fit, wherein said line is parametrically defined according to λ=sθ+ν₀; e) providing said parameters s and ν₀ and said parameters A, B, x₀, y₀, and ε as λ(θ) calibration outputs.
 15. The method of claim 14, further comprising: f) determining a difference e_(i) of said linear angle-wavelength fit at each of said triples according to e_(i)=ν_(i)−(sθ_(i)+ν₀); g) fitting a spline to the set of points (θ_(i), e_(i)) to provide a spline difference fit e(θ); h) providing said spline difference fit e(θ) as part of said λ(θ) calibration outputs. 